Type la supernovae as speed sensors at intermediate redshifts 
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Large scale peculiar velocity (LSPV) is a crucial probe of dark matter, dark energy and gravity at 
cosmological scales. However, its application is severely limited by measurement obstacles. We show 
that fluctuations in type la supernovae (SNe la) fluxes induced by LSPV off'er a promising approach 
to measure LSPV at intermediate redshifts. In the 3D Fourier space, gravitational lensing, the 
dominant systematical error, is well suppressed, localized and can be further corrected effectively. 
Advance in SN observations can further signiflcantly reduce shot noise induced by SN intrinsic 
fluctuations, which is the dominant statistical error. Robust mapping on the motion of the dark 
universe through SNe la is thus feasible to 2 ~ 0.5. 
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Matter distribution of the Universe is being revealed to 
great details by surveys on galaxies, gravitational lensing, 
the thermal Sunyaev Zel'dovich (SZ) effect, CMB, X-ray, 
etc. In contrast, measurements on large scale peculiar 
velocity (LSPV), or bulk flow, are still limited. Measure- 
ments which rely on distance indicators to subtract Hub- 
ble flow [l[ are limited to local universe. Those based on 
anisotropic galaxy clustering in redshift space can be ex- 
tended to cosmological distances i2||. However, this does 
rely on modeling of galaxy redshift distortion, whose ac- 
curacy still requires much improvement in this era ofpre- 
cision cosmology [3]. The kinetic SZ (KSZ) effect [|| of 
clusters is a promising probe Q , although systematic er- 
rors do exist Statistics of the diffuse KSZ background, 
such as the power spectrum jTj and cross correlations 
with other tracers of the large scale structure [8] can be 
measured robustly. However, this background measures 
the peculiar momentum instead and thus probes LSPV 
only indirectly. 

LSPV, as a direct tracer of gravitational potential at 
cosmological scales, is powerful to probe dark matter, 
dark energy Q and gravity In fact, it is indispens- 
able to distinguish between some scenarios on the dark 
sectors We show that it is promising to measure 

LSPV robustly through type la supernovae (SNe la) and 
open a new window into the dark universe. 

Peculiar velocity of a SN at position x shifts the ap- 
parent redshift to 1 + z = {1 + z){l + v ■ x), where 2 is 
the real redshift. It also changes the luminosity distance 
to dL{z) = (ii,(z)(l + 2v- x). Fluctuations induced in its 
flux (with respect to the mean flux at z, instead of at z) 



dUx) =Q{z)v-x; Q{z) = -2(1 



XH 
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Here, H is the Hubble constant at redshift z and x — 
J dz/H is the comoving distance. We have neglected 
the earth motion, whose effect is straightforward to take 
into account. The signal 5"^ ~ v/{cz) is contaminated 
by 5p induced by gravitational lensing^ and S^p^'^°™ in- 
duced by SN la intrinsic fluctuations. 5p|-"dom spatially 
uncorrelated and straightforward to correct. 5p induces 
systematical errors in velocity measurement, however, at 
local universe where z < 0.1, contaminations of Sp are 
negligible {{i6"p)^) > lO"-* > {{Sj^)^)). These properties 
have enabled success in local LSPV measurements fl^ . 

Gravitational lensing induced systematical error in- 
creases with redshift, while the signal decreases. So even- 
tually the lensing induced systematical error will over- 
whelm the signal at intermediate redshifts. For example, 
at 2 ~ 0.5, S^p ~ 0.1% < S^. If not corrected, this 
will prohibit the application of SNe la as cosmic speed 
censors. However, as the main result of this paper, we 
will show that, lensing induced systematical error has a 
intrinsically different pattern to that of 6p and can be ef- 
fectively corrected in the 3D Fourier space. On the other 
hand, shot noise induced by SNe la intrinsic fluctuations, 
which is the dominant statistical error, can be overcome 
by advance in observations. This approach does not rely 
on assumptions on SN la absolute luminosity, other than 
that its rms fluctuation is small, since it only explores 
information imprinted in flux fluctuations. Eventually, 
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^ Dust extinction also causes fluctuations in SN fluxes. This kind 
of fluctuations has opposite sign to that of gravitational lens- 
ing and thus suppresses contaminations in velocity measurement. 
For cosmic gray dust, contaminations can be suppressed by 10- 

50% iia]. 
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high precision mapping of the motion of the dark uni- 
verse to z ^ 0.5 will be realized through SNe la. 

For a narrow redshift bin, fi'"p and i5|; are virtually un- 
correlated, due to the cancellation of positive v and nega- 
tive V and the lensing weighting function. We then have 
{5f{xx)6f{xi)) = A:Wk.{xx,x-2) + Q'^{z)^^(xx.,x-x)- Here, 
the direct observable is the total flux fluctuation (the 
sum of three) 



= SI 
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w^{xi,X2) = {k{xi)k{x2)) and £,y{xi) = {v{xi) ■ Xiv{x2) ■ 
X2) are the correlation functions of k and v{x) -x, respec- 
tively. K is the lensing convergence of a SN la at position 
X and Sp = 2k in the weak lensing regime. Wk{xi, X2) is 
anisotropic with respect to xi —X2 ^- It depends primar- 
ily on 9i2, the angular separation of two lines of sight 
(I.O.S.). The 2D Fourier transform over 6*12 (under the 
flat sky approximation) gives the usual 2D lensing power 
spectrum C^^lxi, 2:2). On the other hand, w^lxi, X2) only 
weakly depends on X1—X2, the radial separation along the 
I.O.S.. Namely lacks small scale power along the l.o.s. 
Correspondingly, the 3D power spectrum P^j^ijt), which 
is the 3D Fourier transform of is highly anisotropic 
in the 3D wavevector k space. Its power should con- 
centrate on a 2D plane of k perpendicular to the l.o.s., 
while its power along the l.o.s. can be significantly sup- 
pressed. On the other hand, P^j-i{k), which is the 3D 
power spectrum of w • has the opposite behavior. Its 
power concentrates along the l.o.s.. This brings hope to 
measure P^jjik) along the l.o.s.. 



THE 3D FLUX FLUCTUATION POWER 
SPECTRUM 

In the appendix, we derive the analytical expression of 
the 3D power spectra P§£,{k) and P^jj{k). These results 
apply to survey volume in a narrow redshift range [z — 
Az/2, z -f A2/2] and sufficiently small sky area flmax x 
^max- We require Az ^ z and 0max/2 <^ 90° and choose 
the coordinate such that the center of the survey volume 
is Xc — (0, 0, Xc). For the small sky coverage adopted, the 
l.o.s. is close to the direction to the survey center, which 
we denote with subscript |j . The directions perpendicular 
to the direction to the survey center are denoted with 
subscript _L. Here we summarize the results. 



^ During preparation of this work, A. Vallinotto et al. |15| and L. 
Hui et al. [lal published works on 3D len sing correlation func- 
tion. Especially in |16| and a pioneer work [TtI by T. Matsubara 
that [T^ pointed out, the anisotropy in the correlation function 
Wk{x\,X2) is explicitly calculated. The result presented in this 
paper is performed in Fourier space, complementary to theirs 
and particularly suitable for LSPV reconstruction that this pa- 
per focuses on. Recently, Hui et al. [Tsj l published their result in 
Fourier space, which is consistent with ours. 
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FIG. 1: The 3D power spectrum variance of SNe la flux fluc- 
tuations, which is the sum of A^(fc) = Q^{z)k^ P^o{k)/{2TT^) 
(solid lines) and Al{k) = Ak^ P3D{k)/{2TT^) (dot lines), at 
z € [0.45,0.55]. k is in unit of h/Mpc. Gravitational lens- 
ing lacks small scale power along the l.o.s., which causes the 
oscillating features. The same reason causes signiflcant sup- 
pression in A|^(fc) at fcy 3> 2tv/Ax, where LSPV can be mea- 
sured accurately in a model independent way. When fcy — > 0, 
A|^(fc) is barely suppressed and A^(fc) can not be measured 
for these modes. A|, begins to increase where k± — > 0, caused 
by its dependence on C^L^^^^. Noticing that A|,(fc) is slightly 
asymmetric with respect to 6 —f n — 6. The ensemble average 
of Aj^ can be calculated from independent lensing surveys. 
After subtraction, gravitational lensing only induces statisti- 
cal error (long dash lines), which is estimated combining all 
^max X Smax patches iu half sky. 



The lensing convergence k is determined by the matter 
over-density 6 projected along the l.o.s. through k{x) = 
f S{xLx)W{xL,x)dxL. The lensing kernel W{xl,x) = 
ifloHQXLi^ — xl/x)/2 where the lens distance xl < x 
and zero otherwise. Wk{xi,X2), Cj^{xi,X2), the 3D mat- 
ter power spectrum Pm{k,z) and the 3D lensing power 
spectrum P^uik) are related by the following equations. 



W{X1,X2) 



(2^ 



.C^(a;l,a;2)e-*'^(*l■^-*^'^^ (3) 



Cr(xi,X2) = / P™( 



and 



I W{xL,Xi)W{xL,X2) , 

— ;zl) 2 "2; I, 

XL xi 



PSoik) = xlAxCi=k^^J{k) 



(4) 
(5) 
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Here, Aa; is the distance interval between z — lS.z/2 and 
z + l^z/2. Ci is Ci{xi,X2) averaged over the source distri- 
bution xi^2- I{k) is an oscillating function oik\\, induced 
by the Fourier transform along the direction of the survey 
center and reflecting the fact that w{xitX2) lacks small 
scale power along the l.o.s.. When fc|| fc±^max and 
g = fc|| A2;/2 > 1, / ~ [sing/gY < 1- Only when fc — > 0, 
/ ~ 1. 

On the other hand, the 3D power spectrum of 6'"p is 



fc4 



(6) 



Here D is the linear density growth factor, f3 = 
din D/dlna and P„i{k,z) is the matter density power 
spectrum at redshift z. The measured flux fluctuation 
power spectrum is the sum of two: 4P;^^ + Q^{z)P^j^. 

We show the two 3D power spectra at z = 0.5 in fig- 
ured] We only show those modes with I = k±Xc > 30, 
where the flat sky approximation holds. The power of 
concentrates along the l.o.s., whereas that of P^j^ 
can be significantly suppressed. Clearly, for modes with 
sufhciently large fc|| , contribution from LSPV outweighs 
that from gravitational lensing by a large factor. The en- 
semble average of P^j^ can be predicted from , which 
will be measured to high precision by upcoming lensing 
surveys on the same sky. By subtracting it from the mea- 
sured flux power spectrum, the systematic error P^jj is 
converted into statistical error and one obtains an unbi- 
ased measure of Pg^ (Fig. [T|). 

One still needs to overcome SNe la intrinsic fluctua- 
tions, which induce the dominant statistical error. 
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Here, fn is the fractional solid angle of k modes used for 
the analysis. 10% rms flux fluctuation (aintr — 0.1) corre- 
sponds to roughly 0.1 mag rms dispersion in magnitude. 
Better calibration of SNe la and improvement over sur- 
vey noise are likely able to reduce a-mtr [3 ■ A factor of 2 
decrease in (Tmtr would relax the requirement of risN by a 
factor of 4 or the sky coverage by a factor of 16. Other- 
wise, successful measurement of LSPV requires (1) large 
sky coverage /sky ^ 1 and (2) high SN number density 
nsN ~ 10~^{h/Mpc)^ . There are no fundamental obsta- 
cles to reach the above goals. The SNe la rate at 2 = 0.5 
is ~ 10~4(/i/Mpc)Vyrl2E So the survey must be 
performed over decades to reach hsn ~ 10~^(/i/Mpc)'^. 
Drift scan using many telescopes of large field of view 
is able to detect all SNe la over a significant fraction of 
the sky. With advance in technology, such survey can 
be feasible and economical . Follow up measurement 
on these ~ 10^ spectroscopic redshifts is also feasible, 
since low resolution spectroscopy with ~ 300 km/s 



suffices. Even for these surveys, the LSPV measurement 
will be still limited to directions with sufficiently large 
however, such measurements contain full information 
on the 3D power spectrum of v (not only v • x), since it is 
isotropic in k space. Furthermore, since the power of the 
velocity power spectrum variance peaks at fc ~ 0.1/i/Mpc 
for the standard cosmology, the bulk information of ve- 
locities can be measured. 

Measuring LSPV through SNe la requires ambitious 
improvement over currently proposed SNe la surveys. It 
thus requires strong justification. Here we briefly address 
some impacts of such measurement on precision cosmol- 
ogy. For cosmology based on general relativity, dark mat- 
ter and smooth dark energy, LSPV measurement helps to 
break parameter degeneracies. For example, the strength 
of weak gravitational lensing increases with dark matter 
density f2m, because of more lens mass, and as, because 
of stronger density fluctuations. Thus there exists a de- 
generacy between r2„j and erg. This degeneracy can be 
broken by the lensing tomography, through which one 
can infer the evolution of the matter density fleld, which 
is more sensitive to flm than to CTg-^ On the other hand, 
LSPV is very useful to break this degeneracy, since the 
strength of LSPV is very sensitive to evolution in the 
matter density field. A simple justification is as fol- 
lows. In the linear regime, from the continuity equation 
(5 4- V • ?; = 0, we have v{k, z) oc fS{k, z) in the Fourier 
space, where / ~ ^rn{z)^'^^ ■ With the measurement of 
d{k,z) (or the power spectrum Pm{k)) from weak lens- 
ing and LSPV measurement, one can directly infer f{z) 
and thus the value of Qm to break the r^m-crg degener- 
acy in weak lensing cosmology. LSPV information is also 
highly complementary to CMB. The strength of CMB is 
sensitive to the shape and amplitude of the initial fluctua- 
tions. LSPV measurement directly provides such kind of 
information and can thus significantly improve the CMB 
cosmological constraints. For other applications, refer to 

i- 

If our universe is indeed fully described by general 
relativity (GR), dark matter and smooth dark energy, 
cosmological probes other than LSPV can already put 
stringent constraints on dark matter density, dark en- 
ergy density and the dark energy equation of state 
In this case, LSPV measurement does not provide funda- 
mentally new information, due to the fixed (and known) 
relations between fluctuations in the metric, density and 
LSPV in such cosmology. However, the dark universe can 
be more complicated. If we are open to the possibilities 
that dark energy may be clustered, with non-negligible 
anisotropic stree, and gravity at cosmological scales may 
deviate from GR, LSPV containes independent and cru- 



^ In the linear regime, the matter density evolution does not de- 
pend on (Tg. But in the nonlinear regime, it does. 
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cial information to that in flucutations of metric and den- 
sity field. LSPV measurement can thus provide valuable 
information, such as the nature of gravity at cosmologi- 
cal scales [lO|. Indeed, to distinguish between some dark 
energy models from some modified gravity models, it is 
indispensable to have LSPV measurements, as shown in 

There are several ways of measuring LSPV. LSPV 
measurement through SNe la is complementary to other 
methods. Among them, a promising one is the kinetic 
Sunyaev Zel'dovich (KSZ) effect of galaxy clusters. On- 
going and planned surveys such as the south pole tele- 
scope (SPT)** plus the dark energy surveys (DES)^ are 
able to measure the LSPV power spectrum to 10% sta- 
tistical accuracy at large scales [5, 6] . Future all sky sur- 
veys can further reduce the statistical error by a factor 
of 3. However, what the cluster KSZ effect directly mea- 
sures is the total KSZ flux cx MgV, instead of v. Here 
Mg is the total gas mass. Some complexities, such as 
feedback, cooling, point source contamination and non- 
isothermal/non-spherical gas distribution, could bias the 
estimation on Mg [§] and thus bias the LSPV measure- 
ment through the cluster KSZ effect. On the other hand, 
LSPV measurement through SNe la is free of such uncer- 
tainties, althrough it may suffer some other systematics. 
For example, if the intrinsic SNe la fluctuations are some- 
how correlated with the large scale structure (although 
quite unlikely), LSPV measurements will be biased by 
such correlation. Nonetheless, SNe la provide an inde- 
pendent, likely clean and potentially powerful method to 
measure LSPV. For all these reasons, surveys capable of 
measuring LSPV through SNe la would be very profitable 
scientifically. 



DISCUSSIONS 

Throughout the paper we adopt the flat sky approx- 
imation, which breaks where I — k±Xc {k± 0). 
However, we will show in a companion paper that, the 
feature that the power of Sp is significantly suppressed 
along the l.o.s. survives beyond the flat sky approxima- 
tion and so as the proposed LSPV measurement method. 

The same technique can be applied to fluctuations in 
the fundamental plane relation of galaxies, for which pe- 
culiar velocity also contributes. Dispersion in the funda- 
mental plane is larger than that of in SNe la. However, 
higher galaxy number density may be reached in shorter 
observation period and the overall shot noise can be re- 
duced. 



As the summary, this paper shows in a proof of con- 
cept study that contaminations induced by gravitational 
lensing in SN peculiar velocity measurements can be ef- 
ficiently and unbiasedly corrected in 3D Fourier space. 
Combining with advance in SNe la surveys, it is feasible 
to measure the large scale peculiar velocity robustly and 
map the motion of the dark universe. 

We thank Neal Dalai, Yipeng Jing and Zheng Zheng 
for useful discussions. This work is supported by the 
one-hundred-talents program of CAS, the NSFC grant 
10533010, 10533030, the NSFC Distinguished Young 
Scholar grant No. 10525314, the CAS grant KJCX3- 
SYW-N2 and the 973 program grant No. 2007CB815401. 



APPENDIX 

We choose the coordinate such that the center of the 
survey volume is at (0,0, Xc). The 3D vector pointing 
to an object is x = xx = {xx±,x\^). We adopt the flat 
sky approximation. Under this approximation, the solid 
angle element dflx ~ d^x± and the volume element dV = 
x^dxdflx ~ x^dxdx±. Also, ~ x, accurate to better 
than 2% accuracy, for 20° x 20° sky. 

The 3D lensing power spectrum P^jj{k) is defined as 



1^ - \ i-r r- - ..dVidV2 , , 
WK,{xi,X2)eyiY>\}k ■ [xi ~ X2\) — — — . (8) 



Plug the Wk.-Ci relation (Eq. [3]) into Eq. [51 integrate 
over (i^xi_L and then d"^!, we have 



x\dxidV2 
V 



Plug into the Limber integral for Cf (Eq. |4]), integrate 
over dxi=i.2 and then integrate over (ix2±, we have 

PSnik) ^ xl/\x [ P^(tl^,ZL)^W^ixL)Iik,XL) 

J XL Xj^ 

= xlAxCi^k^^M) . (9) 

Here we have replaced Pm(%f-) with Pm(%^), since 
Pm varies slowly with k and xi ~ for narrow redshift 
bins. W{xl) = JWx'^dx/ J x^dx is the usual weighted 
lensing kernel. The suppression factor 



I{k,XL) = W-^ixL) 



Ix" W{xL,x) ex-p[ibx]x^ dx 



f " x'^dx 

Jxi 



* http://pole.uchicago.edu/ 

^ http:/ /www. darkenergysurvey.org/ 



Here, b = k\\+ki^ - Xi_ and (■ • • ) is averaged over the solid 

angle d^x±. The term bx in the exponential \sk-x under 
the flat sky approximation. This implies that, an or- 
thogonal Fourier mode is along the actual l.o.s. (instead 
of along the direction to the survey center). The other 
orthogonal Fourier modes are fc^, which is the Fourier 
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transform of angular direction. I(k) is I{k,XL) weighted 
through XL < and its exact definition is given by the 
last relation in Eq. O xi and the distance at 

z — Az/2 and z + Az/2, respectively. 

When XL ^ xi, I{k,XL) varies slowly with respect to 
XL, because W{xltx) in the bracket varies slowly with 
respect to x G and cancels the one in W. When 

xl — * Xu, I increases rapidly and approaches 1. Since 
most lensing contributions come from xl < xi and the 
lensing weighting (W^I) peaks at xl — a;„/2, a con- 
venient approximation is that I{k) ~ I{k,XL = x„/2). 
Numerical evaluation of Eq. [9] shows that this ap- 
proximation degrades for smaller fc, however, even for 
k = (0,0, 0.05)/i/Mpc, it still works to 10% accuracy 
We will adopt this approximation, for its simplicity and 
clear physical meaning. 

At linear scale, peculiar velocities are irrotational. The 
mass conservation equation 5 + V ■ v/a — then tells 
that v{k) = id{k)k/k'^. The 3D (l.o.s.) velocity power 
spectrum is 



[v{xi) ■ X1V{X2) • X2)e 



-ik -(2^1—2^2) 



p ^ dVldV2 



V 



(27r)3 k^ V 
{k^^ + (fc_L • • X2,\_) + 2fc||fc^ ■ 



Pm{k,z).^ 



(fk Pmik)dVidV2 



fc4 II 7 (27r)3 k^ V 



Here, ij_ = {xn^ +X2j^)/2, Aa:j_ = xi± — X2i_ and 
X2,i. — XcAij,. Finally we have 



c? k^ " 



Pm {k, z)kj_ {x^ 
2fc4 



Pm ( ^ 1 z)kikj 



Nonlinearity at small scales could modify the above re- 
sult. For example, shell crossing generates vorticity (ro- 
tational velocity component), whose impact on P^liW 
could be non-negligible where kj_ ^ 0. However, such 
possible modifications do not invalidate the proposed 
LSPV measurement approach. Since the main purpose of 
this paper is to demonstrate the feasibility of measuring 
LSPV through SNe la, instead of developing an advanced 
LSPV model, Eq. [TOl suffices. Furthermore, due to the 
extra factor {kxc)~^, the last term in Eq. [10] is much 
smaller than the second term for relevant k and will be 
neglected elsewhere in this paper. 
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scopes can be constructed with reasonably low cost and 
detect all SNe la over the whole sky. 
P.S. Corasaniti, M. LoVerde, A. Crotts, 
C. Blake. Mon.Not.Roy.Astron.Soc. 
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|http://w ww.astro.ubc.ca/LMT/alpaca| 
http://www.itp. ac.cn/Download/spac2d[)57Arlin-Crotts.pptJ 
For a comprehensive review, refer to the Report of 
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,arXiv:astro-ph/0609591 ) 



